Smearing of phase transition due to surface effect in ferroelectric nanostructures 
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The boundary conditions, customarily used in the Landau-type approach to ferroelectric thin 
films and nanostructures, have to be modified to take into account that a surface of a ferroelectric 
(FE) is a defect of a "field" type. The surface (interface) field is coupled to a normal component 
of polarization and, as a result, the second order phase transitions are generally suppressed and 
anomalies in response are washed out, as observed experimentally. 



Theoretical studies of phase transitions in thin films 
and the corresponding size effects within the Landau the- 
ory [1] have been undertaken since 1950s. Recently, the 
interest to these questions has risen dramatically in view 
of the applications of ferroelectric thin films [2] and a 
discovery of various ferroelectric nanostructures [3] . The 
boundary conditions for thin films were originally dis- 
cussed by Ginzburg and Landau (GL) in 1950 [4] and 
later by Ginzburg and Pitaevskii [5]. It was shown by 
GL that, if the properties of the boundary layer are the 
same as of the bulk, one arrives at the condition that the 
gradient of the order parameter r\ vanishes at the surface, 
Vjjrj = (in zero magnetic field, ft is the normal to the 
surface). Starting from a microscopic theory, de Genncs 
has shown that for a superconductor-metal interface with 
no current and magnetic field a more general boundary 
condition applies, V.377 + r]/S = 0, where 5 is the char- 
acteristic length scale describing the proximity effect [6] . 
These conditions are very general and were introduced 
phenomenologically by Kaganov and Omelyanchouk for 
a surface of a ferromagnet [7] (cf. review in [8]). They 
assumed that the surface energy was oc rj 2 . Kretschmer 
and Binder [9], using the same boundary conditions, have 
taken into account the depolarizing field, which is impor- 
tant for a polarization (magnetization) perpendicular to 
a surface. These boundary conditions are customarily 
used in studies of phase transitions in ferroelectric films 
(see, e.g. [10]). 

The choice of a surface energy dependence on the or- 
der parameter must be governed by symmetry consider- 
ations, as it is for a bulk energy. Surface symmetry is 
different from a symmetry of the bulk, and this is of spe- 
cial importance for ferroelectrics. A surface eliminates 
all symmetry elements that change the direction of the 
normal to the surface. If the order parameter (or one of 
its components) is a vector perpendicular to the surface, 
the surface energy contains a linear "field" term remi- 
niscent of a "local field defect" [11]. The "local field" 
is absent for superconducting or magnetic order parame- 
ters because surface does not break either gauge or mag- 
netic symmetry. A Coulomb dipole field (double layer), 
contributing to the work function [12,13] and structural 
surface relaxation (see, e.g. [14]) are different aspects of 



this phenomenological "field" . If the order parameter is 
not a vector, the "field" term may appear in special cases 
if the surface eliminates appropriate symmetry elements. 
If the field term is absent, some higher order terms, for- 
bidden in the bulk, can be still allowed by the surface 
symmetry. We illustrate the effect of new boundary con- 
ditions for thin film of uniaxial FE. This is relevant also 
to perovskite ferroelectrics that are cubic in the bulk but 
become tetragonal in thin films because of elastic misfit 
with a substrate. Let the surface be perpendicular to 
the polar axis z. The Landau-Ginzburg-Devonshire free 
energy for the bulk between the electrodes 1 and 2 at 
z = is [15,16] 
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where P = P z , and one adds to it the surface energy 
which, generalizing Ref. [7], that can be written in the 
form 
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where a corresponds to a "temperature" -like component 
of the surface energy and w to its "field" component. We 
obtain from Eqs. (1) and (2) the new correct boundary 
conditions for ferroelectrics in the same way as in Ref. 
[7], 
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In CGSE units, one can estimate that a ~ d at , where d a t 
is the characteristic "atomic" length scale, on the order 
of the lattice constant. The electric field at the surface 
[12,13] is on the order of <5> s /d at ~ 1V/A « 10 8 V/cm, 
where q$> s ~ 4eV is the typical workfunction for ferro- 
electrics [2] . The surface field corresponds to a polariza- 
tion, which is on the order of an "atomic" polarization 
Pat = q/dl t ~ 200/^C/cm 2 , so that w ~ P a td a t- We ex- 
pect that the structural relaxation contribution to w is of 
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the same order of magnitude: from first-principles calcu- 
lations at the ideal surfaces of BaTi03 and PbTi03 the 
surface polarization P sur f ~ 10 _1 P at [14]. Usually the 
top and bottom electrode interfaces to the FE film arc 
dissimilar and thickness of the surface polarized layers is 
<~ d at . This is enough to produce a considerable field in 
the film. Indeed, the field in the surface double layer is 
Ei = —4irP sur f <~ — 4nP at . If, for example, we consider a 
short-circuited film with strongly asymmetric double lay- 
ers at the interfaces (so that the voltage drop on the sec- 
ond layer is much smaller than the one on the first layer) 
the surface voltage drop <~ E\d at must be compensated 
by the field in the bulk, Ebuik — —Eid at /l — inP at /l. For 
example, in a film with a thickness I ~ 1000A we find 
Ebuik ~ Patdat/l ~ 300kV/cm, which is a very strong 
field indeed. In a ferroelectric film this field smears sec- 
ond order transitions and may smear a first order phase 
transition. In a free standing FE film the role of elec- 
trodes would be played by a surface conductivity and/or 
atmospheric ions. Note that is a rare case of symmetric 
surfaces the surface dipoles produce no electric field in 
the film and no smearing of the phase transition occurs. 
The polarization P (z) is found from the equation of state 
for (1) and the Poisson equation. Assuming that there is 
no external charge, we have 

AP + BP 3 + CP 5 - gd 2 P/dz 2 = E, (4) 
d (E + AnP) /dz = 0, (5) 

(1/0 Edz = (<pi — ipz) /I = E , (6) 

where Eq is the external electric field. We obtain from 
Eqs. (5) and (6) 

E = E -4n[P(z)-P], (7) 

where the overbar means an average over the film, i.e. 
/ = (1/Z) dzf(z). One probes experimentally not the 
actual surface polarization but the net polarization in the 
sample. Indeed, the experimentally measures the surface 
charge a, which is equal to the normal component of the 
displacement field, a = D z /4n = E^/An + P. Therefore, 
we need to find an equation for the net polarization P, 
which, from the symmetry arguments alone, should have 
the form 

AP + bP 2 + BP 3 + fP A + CP 5 = E + W/l, (8) 

since the terms of all powers in P are allowed in a system 
with asymmetric boundaries. Below, we will solve the 
system (4)-(7) with the boundary conditions (3) and find 
the coefficients in the generic equation of state (8). To 
find the net polarization P, we integrate the equation 
of state (4) over the film with the use of the boundary 
conditions (3): 

AP + BP 3 + CP 5 + (3BP + 10CP 3 )JP 2 



+ (B + 10CP 2 )5P 3 + 5CP5P 4 + C8P 5 
= E Q + R/l, (9) 
R = wi + w 2 — Qi-Pi - a 2 P2, (10) 

where SP = P(z) - P. To find 5P wc subtract Eq.(9) 
from Eq.(4) and obtain 

ASP • K 2 • /v, • A ; • K 5 - gd 2 5P/dz 2 = -R/l, (11) 

where the cumulants K 2 = (3BP + 10CP 3 )(5P 2 - 
IF 1 ), K 3 = (B + IQCP 2 ) (SP 3 - JF^j , K 4 = 

hCP (b~P 4 - JP^j , K 5 = C (5P 5 - SP^^j can be ne- 
glected under conditions used below, A = An + A + 
WP 2 + 5CP 4 w 4tt, since A = (T - T c )/T and for dis- 
placive systems Tq ~ 10 5 K, so we are in a regime where 
\A\ < T c /T < 1. The solution of the equation (11) is 

SP = -R/(Al) + ci cxp [-\(z + 1/2)] 

+ C2 exp[-A(//2-^], (12) 

(~ \ 1/2 -l 
whereA=^/fifJ ~ d at . 

The fact of appearance of the atomic length 1/ A ~ d at 
in our phcnomcnological treatment is important and re- 
quires an explanation. Actually, it has been already 
found in Ref. [9]. Physically, it appears because the 
spatial changes of polarization that we consider lead to 
bound charges and strong electric ( "depolarizing" ) fields 
at the interfaces. As a result, the characteristic length 
of the spatial change of the polarization there is com- 
parable not to the correlation length (the width of the 
domain wall) but to a much smaller length which hap- 
pens to be of an atomic order for typical values of the 
coefficient g. Strictly speaking, the present phenomeno- 
logical derivation is valid for models with ^Jg > d a t ■ Real 
materials are at the boundary of applicability of the the- 
ory (length scale ~ d at ). From the data by Stemmer 
et al. for 90° domain walls in PbTiC>3 [17] we can esti- 
mate 1/A = (g/An) 1 / 2 = 0.4- 0.5A. This suggests that 
PbTi0 3 is not exactly the displacive type system, where 
our estimate applies, 1/A ~ d at ~ 1 — 3 A, but it is compa- 
rable. Therefore, the results below are semiquantitative, 
with the exact values of the renormalized coefficients in 
Eq.(8) to be found from a microscopic theory. 

We are interested in a solution accurate to (leading) 
terms linear in l/l. We find in this approximation 

Cl(2) = ( w l(2) - Oll(2) P) I («1(2) + Xg) , (13) 

and obtain from the condition J SPdz = 

R = * 9 (^f + ^f). d4) 

\ ai + \g a 2 + Xg J 

The averages in Eq.(9) are then easily calculated with 
the use of (12): 5P" = (c™ + c")/(nAZ) for n = 2 - 5. 
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Further simplification is possible if we assume that the 
characteristic length y/g 3> d a t, while a <~ d atl mean- 
ing that we have a small parameter a/Aj « 1. In this 
case the terms -K^2,...,5 contain higher powers of the small 
parameter a/Xg and can indeed be omitted. Expanding 
the results in terms of the small parameter a/Xg <C 1, 



we find 

A = A + a/l, (15) 

W = wi(l - cti/Xg) + w 2 (l - a 2 /Xg), (16) 

a = a 1 +a 2 +(3, (17) 

B = B + AB, (18) 



where (3 = W(w 2 +w 2 )/ (2g 2 X 3 )+5C(wj+w 4 )/ (4 9 4 A 5 ) , 
AB = bC{w\ + w 2 )/(g 2 X 3 l). The terms b and / have no 
effect on a behavior near the phase transition, since they 
are proportional to d a t/l, and can be shown to be small. 
They can be neglected together with the renormalization 
of the coefficient C. The dielectric function is 



A + a/l + 3BP 2 + 5CP 4 ' V ' 

with P determined from Eq.(8) and has a smooth peak 
in all systems with nonsymmetric electrodes. In a sym- 
metric case (W = 0) the second order transition is shifted 
down proportionally to l/l. Indeed, in the Landau theory 
A = (T — T c0 )/T , where T ~ T at ~ 10 5 K in displacive 
ferroelectrics, T c q the phase transition temperature in the 
bulk, and the transition temperature in the film is 

T c - T c0 = -aT /l. (20) 

This result is similar to the one after Kretschmer and 
Binder [9] , but what is new here is that the shift depends 
on surface dipoles via the new term (3 that pushes the 
transition temperature down. Indeed, B is negative but 
small for weak first order transitions as in perovskitcs, 
while C has the usual "atomic" value, C ~ P at . The 
first order phase transition is pushed closer to the second 
order, because of positive renormalization of the coeffi- 
cient B — > B. The shift of temperature of the first or- 
der transition in symmetric case is defined mainly by the 
renormalization of the coefficient A. Indeed, the condi- 
tion of the first order transition A\ = 3B 2 /16C can be 
presented as 

Ai = A Ibulk - a/l + 3BAB/8C, (21) 

where Aibuik corresponds to the transition in the bulk 
and we have neglected the term oc l~ 2 . Let us emphasize 
that if, according to Eq.(21), A\ becomes negative, this 
equation may become inapplicable since the system can 
split into domains. It has been shown in Ref. [15] that 
the ferroelectric transition in a film with perfect metal- 
lic electrodes and a "dead layer" (e.g. vacuum layer) 
proceeds with the domain formation even with dead lay- 
ers down to atomic thicknesses. The imperfect screening 



by metallic electrodes will produce similar results. The 
strongly polarized surface regions that we consider here 
are analogous to those dead layers. The value of A cor- 
responding to the transition with the domain formation 
depends on the materials parameters and can be roughly 
estimated as A ~ —d at / (le 1 / 2 *), where e is the dielectric 
constant in the direction perpendicular to the polar axis 
[15,18]. For perovskites, which are uniaxial due to misfit 
strain, the value of e can be large, so that the condition 
for the domain formation is A < —d at / (^e 1 ^ 2 ) can be 
met in a tiny temperature interval just below the loss of 
stability of the paraelectric phase at the point A = 0. 
Therefore, in experiments with symmetric electrodes one 
would see a phase transition very near A = 0, which is 
close to the transition temperature in the bulk [19]. 

The present theory may, at least qualitatively, explain 
the observed surprisingly strong smearing of the phase 
transitions in thin films, both epitaxial and polycrys- 
tallinc [20,21]. Indeed, we made the estimate of the di- 
electric constants for Bai-^Sr^TiOs (x = 0.3) for thick- 
nesses I = 1200 and 775 nm, which is in good semi- 
quantitative agreement with the recent data [21] (Fig. 1). 
The value of the surface bias field has been in the range 
0.07 — 0.1 CGSE, which corresponds rather accurately to 
our above order-of-magnitude estimates. Slower fall-off 
in the data for the dielectric constant at temperatures 
below the peak may be related to usual presence of do- 
mains in this temperature region. It is worth mentioning 
that there have been various hypotheses put forward to 
explain the observed very strong smearing of the phase 
transition, like a "dead" layer or the relaxor behavior 
[20] . We would like to point out that both of those ex- 
planations are not likely in the case of epitaxially grown 
thin films. In particular, the relaxor model would require 
a presence of strong nanometer scale disorder, whereas in 
e.g. epitaxial PZT films on strontium ruthenate there is 
practically an ideal atomic registry to the substrate [22]. 

The surface polarization discussed above is a special 
case of a polarization due to gradients of a scalar quan- 
tity (concentration c of e.g. oxygen vacancies, density, 
temperature, etc.) and they are accounted for by a term 
like 

f c = - 7 PVc, (22) 

in the free energy, where the coefficient 7 is estimated 
as 7 ~ Patdat [23-25]. This field should be taken into 
account in the case of a film with a compositional pro- 
file (grading) given by e.g. the concentration of one of 
the components of a ferroelectric alloy c — c(z). Gener- 
ally, the other coefficients of the thermodynamic poten- 
tial also become inhomogeneous. There are FE systems 
with intentional concentration profile (graded) [26] and 
unintentional (e.g. with a defect concentration profile) 
that are currently a focus of research. The equation of 
state of the graded ferroelectric film with c = c(z) is 



3 



1000 




200 1 ' ' ' 1 

200 300 400 

T(K) 

FIG. 1. Theoretical dielectric constants for films of 
Bao.sSro.sTiOa with thicknesses I = 1200 and 775nm showing 
complete smearing of the phase transition as a result of the 
surface bias field (see text for discussion). Inset: the data 
from Ref. [21]. 

d 2 P 

A{z)P + BP 3 + CP 5 - g-— - DV\P 

dz z 

dc 

= Ev + Av{P-P)+ 1 —. (23) 
dz 

If it were possible to neglect the inhomogeneity in other 
coefficients of the thermodynamic potential, we have to 
simply add an average value of the "field" jdc/dz ~ 
-fAc/l to the bias field (iSi+iDa)/^ where Ac is the differ- 
ence of concentration through the sample (Ac max = 1). 
Account for the inhomogeneities of the Landau coeffi- 
cients makes the problem more difficult because even 
small inhomogeneities of the coefficients have been shown 
to lead to the domain formation if the bias field is absent 
[18]. As above, the main effect of the surface bias fields is 
smearing out of the phase transition into a monodomain 
state. However, a phase transition with formation of a 
domain structure remains a possibility. For a special case 
of step-wise concentration we have found elsewhere [27] 
that (i) in the case of symmetric boundaries our previous 
results [18] are not changed essentially by the presence 
of the interface field 7 and the soft part splits into do- 
mains. In the case (ii) of asymmetric boundaries the 
domain formation is possible but only for a much larger 
compositional inhomogeneity. 

We thank J.M. Gregg, A. Kholkin, and J.F. Scott for 
stimulating discussions and sharing their data. 
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